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Abstract 



We introduce an effective field theory for the vicinity of a zero temperature 
quantum transition between a metallic spin glass ("spin density glass") and 
a metallic quantum paramagnet. Following a mean field analysis, we per- 
form a perturbative renormalization-group study and find that the critical 
properties are dominated by static disorder-induced fluctuations, and that 
dynamic quantum-mechanical effects are dangerously irrelevant. A Gaussian 
fixed point is stable for a finite range of couplings for spatial dimensionality 
d > 8, but disorder effects always lead to runaway flows to strong coupling 
for d < 8. Scaling hypotheses for a static strong-coupling critical field theory 
are proposed. The non-linear susceptibility has an anomalously weak singu- 
larity at such a critical point. Although motivated by a perturbative study 
of metallic spin glasses, the scaling hypotheses are more general, and could 
apply to other quantum spin glass to paramagnet transitions. 
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I. INTRODUCTION 



Electronic systems with strong randomness and strong interactions have been studied 
in a number of experimental systems including doped semiconductors, metallic alloys, and 
most recently in the doped cuprate and doped heavy-fermion compounds. Some of the most 
interesting physics in these materials arises from the complex interplay of the fermionic, 
charge- carrying excitations and the spin fluctuations. A number of distinct equilibrium 
thermodynamic phases are possible, even at zero temperature (T). In the charge sector, 
we may have metallic and insulating phases (within the insulator we may also distinguish 
further between a Mott insulator, with a true T = charge gap, and a Fermi glass, which 
has localized, gapless, charged excitations). In the spin sector, the ground state can either 
be a spin glass, in which each spin has an infinite-time memory of its spatially random 
moment, or a quantum paramagnet, in which the spin correlations decay to zero in the long- 
time limit. There does not appear to be any fundamental principle constraining the relative 
positions of the transitions in the charge and spin sectors, leading to a rich phenomenology 
of possible T = phases and critical points. 

Previous work has examined the quantum paramagnet phase both in the Mott insula- 
tor PI (where the spin fluctuations can be described by a quantum Heisenberg spin model) 
and the metal 0. Studies of the spin glass phase and its onset have however been mostly 
restricted to the insulating phase. In an infinite-range Heisenberg model of the Mott insula- 
tor an instability of the quantum paramagnet to a possible spin glass was noted [Q. Greater 
progress has been made in elucidating the quantum paramagnet to spin glass transition in 
insulating models of Ising spins in a transverse field (which may be appropriate for insula- 
tors with strong crystalline anisotropy) and quantum rotors [§-0]. The methodology of a 
recently developed Landau theory for this transition in the Ising and rotor models [0] will 
be very useful to us below. 

In this paper, we shall analyze systems in the vicinity of a T = transition between a 
spin glass and a quantum paramagnet occurring while the charge sector is metallic. Our 
motivation for this study comes partly from recent experiments in heavy fermion compounds 
like Yi_rcUxPd3 which appear to show a paramagnet to spin glass transition with increas- 
ing doping (x) in a metallic regime. However, we shall not make comparisons of our theory 
with experiments here, as detailed studies of the vicinity of the quantum transition are not 
yet available. 

We will begin by introducing in Section |T| a quantum field-theory. A, for metallic spin 
glasses; our approach suggests the identification "spin density glass" for such systems. In 
Section |Tl| we will determine the mean field phase diagram of ^ as a function of a quantum 
coupling, temperature, and an external magnetic field. Fluctuations about mean field will 
be studied in Section |V[ first by a perturbative renormalization group (RG) analysis (Sec- 
tion [IV AD , which finds flows to strong-coupling for spatial dimensions d < 8. These results 
will then be used (Section |IV B| ) to motivate a scaling scenario for the strong-coupling region 
in which the critical fixed point involves only static fluctuations induced by the quenched 
randomness. Dynamic, quantum fluctuations are dangerously irrelevant at this static fixed 
point, and their effects are controlled by a crossover exponent —6u < 0. The critical singu- 
larity of the non-linear susceptibility, Xni, is weakened by a positive 9u, thus a non- divergent, 
cusp-like, critical singularity in Xni is possible. This scaling scenario generalizes one pro- 
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posed earlier for insulating Ising and rotor spin glasses [0 which had 9u = 0. Indeed, there 
is no fundamental reason why the insulating spin glasses should not also have > 0. 

Static, or ^ = 0, fixed points have also arisen in studies of some other spin systems. A 
model of quantum rotors in a random field was studied some time ago by Boyanovsky and 
Cardy 0, and their results can be interpreted in terms of such a fixed point. However, 
they did not realize that the crossovers, and positions of phase boundaries, at finite T 
are modified by a positive 9u] the required modification is related to that discussed by 
Weichmann et. al. ||T^ and Millis |[ll|] in a rather different physical context, and will also be 
discussed in this paper. More recently, D. Fisher [1^1 has studied the random Ising model 
in a transverse field in = 1 and shown that the results are consistent with a,h = {] fixed 
point: his scaling results however involve an exponential relationship between energy and 
length scales, unlike the more usual power law relationship which we shall find. Finally, 
in very recent work, Kirkpatrick and Belitz |jl3| have proposed a scaling scenario for the 
metal-insulator transition which appears to have many similarities to our results below on 
the metallic spin glass to paramagnet transition. 



II. EFFECTIVE ACTION 



An initial analysis of metallic spin glasses was performed some time ago by Hertz ||14|| , 
although he did not focus on the vicinity of the T = quantum transition. We will study 
models, similar to those in Ref |T^ , described by the following class of Hamiltonians: 

H ^ ^ tijCiiy,Cia ^ ^ J ijSifjiS j -\- Hint, (^■■^) 

i<j,ct i<jtl^ 

where Qq, annihilates an electron on site i with spin a =T, i, and the spin operator Si^ = 
4a'^a/3Cj/3/2 with 0"'^ the Pauli matrices. The sites i,j lie on a (i-dimensional lattice, the 
hopping matrix elements tij are short-ranged and possibly random, and the Jj^j are Gaussian 
random exchange interactions, possibly with spin anisotropies. The remainder Hint includes 
other possible short-range interactions between the electrons: we constrain them so that the 
ground state of H is metallic. A version of H with infinite-range hopping and exchange was 
studied recently []T3| using a static ansatz for the order parameter but with no additional 
approximations. We will provide below a theory which includes dynamic, quantum effects 
and also applies to models with finite-range interactions in finite dimensions. Our analysis 
of H will be similar in spirit to the Stoner model approaches to the appearance of spin- 
density- wave order in clean metallic systems , except that we now consider condensation 



into a density wave with a random orientation of spins, or a "spin density glass" . 

We now derive a low-energy field-theory for H in the vicinity of the spin glass to para- 
magnet transition. The procedure is similar to that of Ref 0. The metallic nature of the 
system expresses itself mainly through the modification of a single term, which, however, 
has important consequences. We begin by defining the spin glass order parameter 

Qf,{x,Tr,r,)= Y: Sl{n)Sl{T,), (2.2) 
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where a,b = l...n are replica indices (the hmit n — is taken at the end), ti,T2 are 
Matsubara times, and is a coarse-graining region around x. At T = 0, the Edwards- 

Anderson spin glass order parameter, qea is the expectation of the replica-diagonal part of 
Q in the limit \ti — T2I — * 00; however, it is necessary to retain the time dependence and all 
replica components of Q to capture the full structure of the field theory The effective 
action for Q is obtained by averaging over the Jij randomness in H after introducing replicas, 
decoupling the resulting 8 fermi term by a Hubbard Stratonovich field Q, and integrating 
out the fermions. This procedure has been carried out in Ref for the quantum rotor spin 
glass, and in Ref |TB] for spin density wave formation in metallic systems. It is simple to 



combine these methods and we omit all intermediate steps. The final effective action. A, is 
expressed in terms of a shifted field Q ^ Q ~ C^y5"''^5{Ti — T2) where the subtraction only 
modifies the uninteresting short time behavior, and the constants Cf^i, are chosen so that 
the resulting Q is small at low frequencies near the critical point 0]: 

A-U^U^ri: r,Q-(., r,r)-±J ...... £ 

+ - / dTidT2 \yQtu{x,Ti,T2) - - / dTidT2dT3 Y Qfui^, n, T2)Q%{x, T2, T3)Q'pl{x, T^, n) 

+1 [drYb Q7^i^.r,T)Q;,l{x,T,T) + v Q^;^{x,t,t)QZ{x,t,t)]] 

-^Jd'^J Y Q';;i^,n,n)Q''Z{x,T2,T2) + ■■■ . (2.3) 

ab^u 

We have only displayed the small subset of terms which will be important near the critical 
point. We have allowed a /i dependence in to reflect possible spin anisotropies; the less- 
important dependence of other couplings has been suppressed. The metallic nature of 
the system is reflected in the second term which has a long-range l/(ri — ra)^ interaction 
in time; the power-law decay is a consequence of the gapless particle-hole spin excitations 
which lead to the dependence 

m^in)SUr2))] ~ V(n - r2)-^ (2.4) 

in any metallic paramagnet |l|,[l6| (the angular brackets represent an average over quantum 
and thermal fluctuations, and the square brackets are an average over randomness). This 
behavior is of course only valid at large |ri — r2 1 and we cut it off at short time differences such 
that its Fourier transform is ~ — ji^l for small ut. The remaining terms in A are identical to 
those obtained in Ref for the rotor model. These are the most general terms, of low order 
in Q, which are local in space-time and consistent with underlying symmetries. In particular, 
the time arguments of Q associated with different replica indices must always be integrated 
independently because the disorder is static. "Quantum-mechanical" interactions occur only 
within the same replica, allowing a gradient expansion about the equal-time point for such 
terms. A more detailed discussion of these criteria can be found in Ref 0] . The particle- hole 
continuum will induce non-local corrections to these terms, but none are as important as 

that in the term linear in Q. The symmetries also allow a 'mass' term ~ ^(^^^(x, ri, ra)^ , 
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but such a term is redundant as it can be removed by the shift Q ^ Q — C^y5°''^5{Ti — T2); 
the shift has a delta function in time and thus does not modify the long-time, low frequency 
behavior that we are interested in. 

We now discuss the physical significance of the couplings in A; for more details the 
reader is referred again to Ref ^ and we highlight only the main points here. First note 
that there are more terms than coupling constants, but it is easy to check that rescalings of 
X and T always allow one to reach the form chosen. The coupling multiplies what turns 
out to be "thermal operator" which tunes the system across the transition. The important 
non-linearity is the cubic k coupling which is induced by disorder effects and involves no 
exchange of energy between the Q fields; there is a I/k in the linear term to ensure that the 
bare Q propagator is independent of k. The only terms involving energy exchange are the 
quadratic u and v terms which represent the "quantum-mechanical" interactions between 
the fermions. Finally, the term in A represents disorder fluctuation effects and arises 
from fluctuations in the local position of the critical point as determined by r^. 

III. MEAN FIELD THEORY 

We now consider the mean field (or tree-level) properties of A. We will only consider 
two extreme limits of the spin anisotropy — Ising-like, when ri ^ r2,r3, and Heisenberg-like, 
when ri = r2 = r^. We will drop the vector fi index except where necessary, and represent 
the effective number of components by M (= 1 for the Ising case, and = 3 for the Heisenberg 
case) . 

In Ref we found that, as in the classical spin glass, a theory with only a cubic non- 
linearity was replica symmetric, even in the spin glass phase; replica symmetry breaking 
appeared only upon including a certain quartic coupling which was formally irrelevant at 
the quantum critical point. Much of this structure carries over unchanged to the metallic 
spin glass, and we shall not dwell on it here. We will restrict our considerations to the 
replica symmetric theory which contains all of the dynamic effects which lead to the im- 
portant differences between metallic and insulating spin glasses. We therefore make the 
following replica-symmetric, x- independent ansatz for the mean field value of Q: in Mat- 
subara frequencies (which are integral multiples of 271 kBT/h as usual) we write 

Q^'iujuuj^) = /3'5.„o5c.„og + /35"'5..i+.„ox(^^i), (3.1) 

where j3 = h/kBT, and x is the local dynamic spin susceptibility (we will absorb a factor 
of kB/h into T from here on). The first term is the spin glass order parameter — note 
that it is independent of replica indices, and, therefore, the replica diagonal and off-diagonal 
components of q are equal and q = qEA- In a theory with only a cubic non-linearity (to which 
we shall restrict ourselves here) the equality of all the replica components of q holds at all 
T; upon including higher-order terms in A, the equality between all the replica components 
persists at T = 0, but there are thermal corrections at any non-zero T which distinguish the 
diagonal and off-diagonal components and which also break replica symmetry 0. 

We now insert ( |3.1| ) into A and determine the saddle point. This determines the spin 
glass to paramagnet phase boundary at r = rc(T) and the order parameter 
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^ 1 (re(T) - r)/{K(u + Mv) for r < r,(T) 2) 
^1 otherwise ' 

with 

r^(T) = r, - c(m + Mv)T^/^, (3.3) 



where c = y7r/2^(3/2), rc = rc(0) ~ Af,/^ is dependent on the large frequency cutoff (see 
Fig p. The local, dynamic spin susceptibility has imaginary part 

^"(cu) = --ImV-tco + A = (3.4) 

Note that the spin fluctuations are gapless, and the crossover from paramagnetic to crit- 
ical fluctuations occurs at a frequency scale A; A also determines the correlation length 
^ ~ A~^/^ that appears in spatial correlation functions, which can be found as Gaussian 
fluctuations around the saddle point. The value of A has a rather complicated dependence 
on T and r and we describe its limiting behavior in the five different regimes of Fig ^] — 
there are smooth crossovers between these regimes. Within the spin glass phase, we have 
A = 0, x" ~ sgn(co')a;^/2 everywhere and there are no crossovers in the present approxima- 
tion. However, we expect that there will be a crossover between a region characterized by 
the quantum ground state (I) to a region dominated by thermal, critical fluctuations (II) 
and such a crossover boundary has been shown in Fig |I|. In the paramagnetic phase the 
present approximation is much richer, and shows all the expected crossovers. The scale A 
is determined by the equation 

A = r- (n + Mt;)^5](|cj| + A)i/l (3.5) 

Solution of (|3.5|) yields four different regimes of behavior (II- V) which we describe in turn. 

(II) |r — rc(T)| <^ [u^Mv)^T'^: this is closest to the phase boundary, and is the region with 
classical fluctuations. We have A = ((r — rdT)) / {T{u + Mv))Y — note that A depends on 
the square of the distance r — rc(T) from the transition (as shown in Ref 0, this is crucial 
for obtaining the classical exponent v = 1/2). 

(III) (|r — rc\/{u + Mv)Y^^ <C T: this is the 'quantum-critical' region in that T is the most 
significant energy scale and the system behaves as if it is at the critical coupling r = Tc- 
Here we find, to lowest order in u^v that A = c{u + Mv)T'^^'^. Note that the expected 



quantum-critical scaling A ~ T is violated. This is a consequence of the fact that all 
T dependent corrections are controlled by quantum interactions m, v which are irrelevant at 
the critical point (see below) — in other words > spoils the naive A ~ T scaling. A 
similar interpretation can be given to the position of the phase boundary rc(T). 
(IV) r — ^ T ^ ((r — rc)/(u + Mv)Y^^ and (V) T <^ r — r^: these are the 'quantum- 
disordered' regions in that T dependent corrections are secondary and to leading order in 
T, M, V we have A = r — r^. The subleading terms in A are different in the two regimes: 
in regime IV, A(T) - A(0) = c{u + Mv)T^/'^, while in regime V, A(T) - A(0) = {u + 
Mv)nT'^/ (Gi/r — fc). This subdivision of the quantum-disordered region is similar to that 
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found in a different context in Ref |TT[ , and is also a consequence of tlie dangerous irrelevancy 
of u, V. 

All of the above crossover boundaries and exponents are of course characteristics of the 
mean field theory, which can, in general, be modified by fluctuations — we will indicate the 
nature of these modifications in the discussion in Section [IV B[ 



A. Phases in a magnetic field 

We complete our discussion of mean field theory by discussing the effect of an external 
magnetic field, if^ on A. The additional terms induced by H can be determined following 



Refs ||T^ which examined the effects of H on spin-density wave formation in clean systems- 



by this method we found 

ab 

--yd xdT}_^ \ iaie^^xHx — 



+ a2HxHxQZ{x,T,T) 



Tl=T2=T 



+ a,{H,H, - HxHx6^,)Q';Zi^, r, r) j . (3.6) 

The field has several different competing effects. The first term, proportional to the coupling 
g, is the static paramagnetic susceptibility of the fermions which polarizes the spins along 
the field, and which always dominates at small H. The a terms account for the precession 
of the spins in the plane perpendicular to H and the energetic contribution of quantum 
fluctuations about the static spin directions: as in clean antiferromagnets [|IU| we expect 



these terms to prefer magnetic order in a plane perpendicular to H (and so as > 0). We 
now consider some cases separately: 

(z) Ising spins, H along easy-axis — only the term proportional to g in ( p.6| ) need be consid- 
ered as the a terms are never important. The finite field phase diagram, the field dependence 
of observables, and the position of the Almeida-Thouless boundary |^ , are essentially iden- 
tical to that for the insulating Ising model considered earlier 0, and will therefore not be 
considered here. The only difference in the metallic case is that no logarithms are present — 
e.g. the free energy at r = rc and T = that depended on H as H^^^ / \og^^^ Hin Ref 0, 
here varies as 

(a) Ising spins, H perpendicular to easy-axis — now the g term in ( |3.6| ) couples only to non- 
critical components of Q and is not important; after integrating out the non-critical Q we 
find that the main consequence of the a terms is to induce a shift ri — > ri + a'H'^ in the 
position of the critical point. 

[Hi) Heisenberg spins — in finite field we now have to allow for the possibility of spin glass 
order appearing in the plane perpendicular to H; the onset of this order is the Gabay- 
Toulouse transition at H = Hqt 0- Let us take a field H pointing along the /i = 3 
direction. The subsequent mean field theory is most convenient in a circularly-polarized 
basis for the vector components of Q: we take Q33 = Ql, Qu = Q22 = {Q+- + Q-+)/2, 
Q12 = —Q21 = i{Q+- — Q-+)/2, and all other vector components of Q = 0. We make 
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the same ansatz as in ([Ol) for the frequency and rephca dependence oi Ql, and Q-+ 
by introducing the quantities ql, xl, 1+- etc.. It is then not difficult to solve the result- 
ing mean field equations. It is slightly more convenient to approach the Gabay- Toulouse 
boundary from the Gabay- Toulouse phase with spin glass order in the transverse direction: 
g+_ = q ^. ^ 0. The solution of the mean field equations for this case are 

= = It 

= = (l^^l + iaxRu)^!'^ 



K 



Ql 



^^(,^) = _1(|^| + A)V2, (3.7) 



K 



where the Gabay Toulouse spin glass order parameter qx and the frequency scale A are 
determined by the solutions of the two equations 

A = r + a^H' + (u + v) - ^ Y.{\uj\ + A) + 2v (yKqr - ^ ^d^l + ^o^lHuf'^^ 

+ (u + 2v) (nqr - Ed^l + ioiiHuf'^^ . (3.8) 



= r + (^2 - a5)H^ + v 



The Gabay- Toulouse boundary is determined by imposing the condition = on these 
two equations, which gives us a line H = HcTif) in the r — H plane. The result of such a 
computation at T = 0, is shown in Fig ^ For small H, the first term in (|3.6|) dominates 
and we find Hgt ~ {"^c ~ rY^^- For large if, the a terms take over, and for 03 > we find 
that Hgt turns over and extends to r > as Hgt ^ {r — Vf)^!'^ (see Fig for sufficiently 
negative 03 this turn over will not occur. 



IV. FLUCTUATIONS 

We begin by a perturbative RG analysis of fiuctuations, which will unfortunately not be 
of much direct utility as there is runaway fiow to strong coupling below c? = 8. Nevertheless, 
the structure of this analysis will help motivate a general scaling scenario which we will 
describe subsequently. 



A. Perturbative RG 

The perturbative RG analysis is quite similar to that of Ref 0. The main difference at 
tree level will be that the dynamic exponent z is z = A rather than z = 2. This difference has 
the important consequence of now making the u,v couplings dangerously irrelevant, which 
in turn leads to a positive 9u- 
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The RG begins with the rescahngs 



x' = xls , r' = r/s^ , i! = ts-' , = gs^^^^+^^-^+^^Z^ (4.1) 

The exponents ri have their usual meaning, while Q is introduced to allow for violations 
of hyperscaling: we will have 6* > causing t to flow to 0, and behave as a dangerously 
irrelevant variable. The irrelevant coupling t and its exponent ^, should not be confused 
with the couplings controlling quantum-mechanical effects and the exponent which will 
be discussed momentarily. At tree level (or equivalently, at the gaussian fixed point), the 
above rescalings leave A invariant provided we modify the couplings 

u' = us^-'-'^ v' = vs^-'-"^. (4.2) 

and choose the exponents 

z = 4 ri = 9 = 2. (4.3) 

Thus the cubic non-linearity k becomes relevant for d below 8, and the rescaling of r gives 
us the gaussian exponent u = 1/4; note the correlation length is given by ^ ~ A-i/^ in the 
notation of Section |T|. The most important point, and the key difference from Ref 0], is 
that u and v are now irrelevant with exponent —2. As these are the only couplings associated 
with quantum effects, we introduce a new crossover exponent, —6u which will control the 
'dangerously irrelevant' consequences of quantum fluctuations. At tree level we clearly have 
6u = 2. Finally, note that for d above 8, n also becomes dangerously irrelevant about the 
Gaussian fixed point. 

It is straightforward to extend the above analysis to include one-loop diagrams. Using 
the diagrams discussed in Ref 0], we find (s = e^) 

z{i)=A + 8K\e) , r] = 2K\i) , 9 = 2, (4.4) 

and the flow equations 

^ = zrii)-an\i) ; = '-^-i^) + ^-'iO- (4-5) 

We have absorbed various phase-space factors into the couplings (see 0), and a is an 
uninteresting positive constant. The irrelevant couplings u, v were set equal to at the 
outset. There is no stable fixed point of ( [4.5|) for real n below d = 8. Above d = 8, the 
Gaussian fixed point is stable, but its domain of attraction is limited to a region which 
vanishes as d approaches 8 from above. For d < 8, and for all physical initial conditions, 
the coupling k flows to strong coupling, making quantitative computation of exponents 
impossible in the present approach. 



B. Scaling Hypotheses 

We will now discuss a non-perturbative scaling scenario for quantum spin glasses, assum- 
ing that the structure of the dangerously irrelevant variables remains similar to that found 
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in the perturbative analysis above. We will consider a static strong-coupling critical theory 
with two dangerously irrelevant directions: one associated with a coupling analogous to t 
which controls disorder fluctuation effects and has exponent —6, and a second associated 
with dynamic, quantum mechanical effects (couplings u,v) and exponent —9u- In the pre- 
vious analysis of insulating spin glasses ^ only t, the first of these dangerously irrelevant 
couplings, was present; our present scaling relations reduce to the earlier ones upon putting 
6u = 0. As we shall see below, a positive 6u has new and important physical consequences. 
Although the analysis below is clearly motivated by our mean field theory of metallic spin 
glasses above, there is no fundamental reason why the insulating models considered in Ref 
should not also have 6'„ > 0. 

It is helpful to discuss non-perturbative effects by considering the scaling behavior of 
observable correlation functions. Among two-point correlators of Q, there are three inde- 
pendent observables 0: in the paramagnetic phase these are the spin glass susceptibility, 
G, 

G{X -y,Ti- T2, Tg - T4) = [{Si^^in) S j f,{T2)) {Si^ir^) S j^ir^))] 

flU 

= H^j;hT) ^ ((«(^'-i'-3)gi(l/,r„r4))) (4.6) 

(the double angular brackets represent averages with a replicated, translationally invariant 
action like A)\ the quantum mechanically disconnected correlator, G"^, 

G'^{i - j, Ti - T2, T3 - T4) = [{Si^{Ti)Si^{T2)) {Sju{T3)Sju{T4^))] - subtractious 

= hm — ^ Y: ( ^1' ^2)Qt{y, r-s, u))) (4.7) 

and the connected correlation function, G"^, 

Gl^ypA^ - j, Ti - ta, T2 - n, T3 - n) = [{Si^{Ti)SiuiT2)Sjp{T3)Sja{T4))] - subtractious 

= li^f?, - E {{QTui^^ r2)Q;:{y, ra, r4))) - • • ■ . (4.8) 

The non-linear susceptibility, Xni, is given by the integral over space and time of G^. The 
correlator is non-zero only because of the "quantum interactions" u, v , and therefore 
carries a prefactor of u, v; in contrast G and G'^ are non-zero even in a purely static theory. 
Under the rescalings ( [4.1| ) we may conclude from arguments similar to those in Ref ||^ that 
G'^ and G scale as 

G'^(x, r, r) ~ (4.9) 
^(x, r, r) ~ , (4.10) 

for fixed t/x^ at criticality. The scaling dimensions of G and G'^ differ because G carries a 
prefactor of the dangerously irrelevant variable t, while G'^ does not. Finally, G'^ carries a 
prefactor of t, and an additional prefactor of the irrelevant quantum interactions: hence 

G'{x, r, r, r) ~ x~^d+2z+e^-2+v) _ 
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Indeed, one can consider the 3 equations ([4.9|) - ([4.11|) as the definition of the three indepen- 
dent exponents rj, 6, and 6u- In the previous analysis [0], 6^ = and hence G and had 
the same scahng dimension. By taking the spacetime integral of ( |4.8| ), we can deduce that 
the non-linear susceptibility Xni behaves as 



Xnl 



-{2-r)-eu)v 



(4.12) 



near the T = quantum critical point. For sufficiently large Xni need not diverge. 

The existence of a static critical theory, and the associated positivity of 9^ , has important 
consequences for the finite T behavior away from the critical point. Recall that T only 
appears as a finite-size length, 1/T, along the time direction, and hence scaling ^ implies 
that its scaling dimension is z. However some power of a combination of interactions like 
M, V must appear in any frequency scale and hence we expect that the naive scaling of finite 
uj correlators as functions of uj/T will now be modified. A related modification of naive 



scaling has been discussed in Refs [0,|TT|] for some clean systems, and we will now present 
a similar analysis. It is useful to consider a simple model of the renormalization group 
fiows near the quantum-critical point at low temperatures. Let us move away from the 
quantum critical point (r = r^, T = 0) by perturbing the system along the single, relevant 



eigendirection by the amount r 



and along the least irrelevant eigendirection which 



involves terms with frequency exchange by the amount u. For small r, u, and T we expect 
fiow equations like 



dT{i) 

di 
dr{i) 

di 
du{i) 

de 



zT{i) 

^irii)-r,)+uf{Tii)) 



(4.13) 



where /(T) (/(O) = 0) is some function arising from thermal occupation of the short distance 
modes of the order parameter fluctuations which are begin integrated out. The key property 
of ( [4.13|) is that a T dependence is induced into the flow of the relevant coupling r only via 
the irrelevant coupling u. The integral of ( |4.13| ) is 



r{i) - r, = (r - r,)e'/'' + ue'/" [ dee-^'^+'Z"^'' f{Te''') 



(4.14) 



(as is customary, we have abbreviated r(0) = r, u{0) = u and T(0) = T) . We now change 
integration variables io Q = Te^^ , and integrate to the correlation length ^ = e^^^ at which 
r(£*) — Tc = 1 to obtain 



z Jt C 



/(C) 



(4.15) 



It is now possible to deduce scaling properties provided it is permissible in the critical region 
to set the lower and upper limits of the integral in ( [4.15|) to zero and inflnity respectively. 
As /(T) represents thermal contribution of short distance modes we expect it to vanish as 
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T — s> 0; these modes however do mix with the particle-hole continuum of the metal, leading 
us to expect a linear density of states at low energies even at short distances, and therefore 
/(T) ~ for small T. In the opposite large T limit, all modes must become classical, and 
therefore f{T) ~ T. For these asymptotic behaviors in /(T), the limits on the integration 
can be extended provided zu < 1 + duV < Izv. We then obtain at r = Tc but T finite 
urp{i+e^y)/z^ In this same region, a similar reasoning implies that the local dynamic 
spin susceptibility will scale as 

x''(c.)=.'-»--'V..^(^_|i_). (4,16) 

for some universal scaling function 0, and non-universal constant K. At tree-level this gives 
us a frequency scale ~ T'^/^ which is consistent with the results of Section |T|. Similarly, the 
position of the finite temperature spin glass to paramagnet boundary (at r = Tc(T)) will 
scale near the quantum critical point at r = rc and T = as 

r, - re(T) ~ mT(^+^"'')/"'^; (4.17) 



Again, this agrees with the tree- level result ( p.3|) . A very similar result applies to the 
boundary between regions III and IV of Fig |l] which occurs at r - rc ~ uT'-^+^^''^^'^ , while 
the boundary between regions IV and V is at r — ~ j^i/zf. 

All of the results discussed so far in this section have been obtained using only rather 
general scaling ideas. In particular, they do not rely on the particular form of the action A. 
We will now obtain a few results which do rely on explicit features of A, and their validity 
is therefore somewhat more questionable. 

A simple argument can be given to fix the value of z, using the manner in which time 
dependence enters into A. Consider a correlator of the Q fields in which all external fre- 
quencies have been fixed at the same frequency u. As the critical field theory is static, and 
because the Q field is bilocal in time, uj will act simply as an external source which shifts 
the value of the "thermal" coupling r — > r + as is apparent from the first two terms 
in A; for insulating Ising and rotor models the corresponding shift is r — > r + cj^. As the 
scaling dimension of r is l/i/, this gives us the scaling relation 

_ f 1 metallic spin glasses ,^ 

1 1 /2 insulating Ising and rotor spin glasses 

We emphasize that both results rely on the assumption of a static critical theory; this 
assumption was not made in the analysis of Ref . We also note in passing that the present 
argument fixing the value of z cannot be applied to the random-field quantum rotor model 
of Ref (which also had a static critical point), because the same external frequency a; 
does not flow through all internal propagators in this case, and some propagators are always 
at zero frequency. 

We now ask whether there is a classical statistical mechanics field theory which is also 
described by the static critical point postulated above. We are only able to answer this 
question within the confines of perturbation theory: a perturbative expansion in k, suggests 
that the relevant field theory is that describing singularities along the imaginary field, ih, axis 
in a (i-dimensional randomly diluted Ising ferromagnet . The latter model has a Yang-Lee 
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edge singularity ||22| at the same value, h = h'^ as the non-diluted Ising ferromagnet p3| . 
Note however that = in random Ising ferromagnets with an unbounded probability 
distribution for the local randomness. It has been argued that there is critical field, 



h = he, such that for < h < he, the zeros of the partition function are analogous to the 



localized states in the band tail in Anderson localization. (The 'Griffiths effects' leading 
to this region also have a parallel in the paramagnetic phase of the quantum spin glass.) 
The singularity ai h = he is then analogous to a mobility edge |^ 
h = he > 0, called the 'pseudo Yang-Lee edge' in Ref pl| 
theory for this singularity is |5^,2^J2T|: 



It is this singularity at 
that interests us here. The field 



A 



YL 





+—y 



(4.19) 



where 0" is the replicated order parameter for the Ising model. This field theoretic model 



was also considered earlier by Parisi and Sourlas |2^, who argued that for k imaginary, 
Ayl describes the statistics of lattice animals. As we will argue shortly, the perturbative 
RG equations for Ayl are given precisely by ( [4.5|) , and a perturbative fixed point with k, 
imaginary can indeed be obtained in the 8 — d expansion. However, in this paper we are 
only interested in the case of k real, which also describes the 'pseudo Yang-Lee edge' |^ in 
a random Ising ferromagnet. 

Now we discuss the perturbative connection between models defined by Ayl and A. 
Consider the Feynman graph expansions with the action Ayl for the correlators 

1 



GYLix - y) 



- y) 



n 



a 

1 



y) 



n[n 



rix)<l>\y) 



(4.20) 



Compare this with the FejTiman graph expansion with the action A of zero frequency cor- 
relators G and G"^ respectively. It is not difficult to show, term by term, that these two 
expansions are identical to all orders in k and to leading order in t. The fact that A involves 
a matrix field (two replica indices) while Ayl has scalar field (one replica index) does not 
affect any of the multiplicity factors associated with any graph; to leading order in t, all 
relevant graphs were tree graphs before averging over the disorder in that corresponds to 
the 1/t^ vertex, and none of these graphs have numerical factors associated with summation 
over replica or vector indices. The equality of these perturbative expansions suggests, but 
does not establish, that the perturbatively inaccessible static fixed points of A and Ayl may 
also be identical: if so, any scaling relations satisfied by Ayl should apply also to A. 

The non-random Yang-Lee edge problem has a simple scaling structure [0] — there is a 
scaling relation between the exponents rj and u as the order parameter (p is also the "thermal 
operator" . The simplest scaling hypothesis for the random case is that the identification of 
(f) as the thermal operator continues to hold. This gives us the scaling relation 



1 d-e + 2-r] 



(4.21) 



Numerical tests of this scaling relation in the randomly diluted Ising model would be quite 
useful (numerical studies of the full quantum spin glass problem are expected to be much 
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more difficult). When combined with ( [4.1^ ), ( [4.21| ) leads to a scaling relation between 9, 
and 7] which is very similar (or identical if zv = 1) to one considered recently by Kirkpatrick 
and Belitz for the metal insulator transition. 



V. CONCLUSIONS 

This paper has proposed a quantum field theory, defined by the action A, for the low 
energy properties of metallic spin glasses in the vicinity of a T = transition between a 
metallic paramagnet and a metallic spin glass. The mean field phase diagram of the model 
as a function of a quantum coupling, temperature and applied magnetic field was described. 
The phase transitions and crossovers in this phase diagram were argued to be characteristic 
of a zero temperature, static critical theory containing no dynamic quantum fluctuations. 
Quantum effects were shown to be dangerously irrelevant and controlled by a crossover 
exponent —6^ = —2. 

Next an attempt was made to extend these results beyond mean field theory, but found 
runaway flows to strong coupling for all spatial dimensions below d = 8. Nevertheless we 
used the insight gained from the mean field theory to propose a set of scaling hypotheses. 
We assumed that the true critical theory also contained only static, randomness induced 
fluctuations, and the exponent 6^ controlling quantum effects took an unknown positive 
value. This had some important observable consequences, similar to those found in the 
mean field theory: 

(z) The non-linear susceptibility had a weaker singularity at the critical point than might 
have been suggested by the usual scaling arguments. In particular, for a large enough 9^ a 
non-divergent cusp-like singularity was also possible. 

(m) In a simple model of the renormalization group flow equations, finite T dynamic response 
functions scaled as functions of uo/T^^^^^ , rather than the usual scaling as functions oi uj/T. 
(Hi) Exponents associated with various crossovers in the vicinity of the T = critical point 
were modified by 9u- 

While these results were directly motivated by our analysis of metallic spin glasses, it is 
possible that some of the scaling ideas are more general and could apply also to insulating 
Ising and rotor spin glasses and other T = transitions in random quantum systems. 
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Note added: In a recent preprint [27|, Sengupta and Georges have considered a model 



closely related to ours, and obtained results in general agreement with the mean-field theory 



of Section III 
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FIGURES 



FIG. 1. Phase diagram of the action A (Eqn ( |2.3D ) as a function of temperature T and r which 
measures the strength of quantum fluctuations. The full line is the only phase transition and 
dashed lines denote crossovers between different regimes, which are described in the text. 

FIG. 2. T = phase diagram of A for the Heisenberg case, in a field H. Hqt is the 
Gabay- Toulouse boundary, qi = [(5^)^] for ^ along the field direction, and similarly qx = [(S*^)^] 
for II perpendicular to the field. 
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Sachdev, Read, and Oppermann Fig 1 




Sachdev, Read, and Oppermann Fig 2 



